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I describe the current status of the theory of tttt scattering, reviewing in particu- 
lar recent work on the numerical solution of Roy equations and on the matching 
between these and the chiral representation. I discuss numerical results on the 
scattering lengths and other threshold parameters. 



1 Introduction 

The study of tttt scattering is a classical subject in the field of strong interac- 
tions: it is a scattering process that may occur at very low energy, far away 
from the resonance region, and it only involves the quasi-Goldstone bosons 
of the chiral symmetry of strong interactions. If we stay at low energy, we 
could expect its dwnamics to be strongly constrained by chiral symmetry: as 
Weinberg showed 111, in the chiral limit, the scattering amplitude has to vanish 
when the momenta of the pions tend to zero. In the real world quarks are 
not massless, and consequently also the would-be Goldstone bosons acquire 
a small mass. The chiral argument given above implies that the amplitude at 
threshold must be proportional to the square of the mass of the pions EJ: 

where stands for the scattering length in the isospin / channel with angular 
momentum l. 

The formulae above are valid only at leading order in a series expansion in 
powers of the quark masses: the nexL-to-leading order corrections have been 
calculated by Gasser and Leutwyler and even the next-to-next -to-leading 
order corrections are now known u. The latter had been calculated already 
before the previous edition of the "Chiral Dynamics" series of Workshops, and 
the interested reader may find a report of the situation at that timeJn the 
contribution by Gerhard Ecker to the previous Workshop Proceedings □. The 
main difference between that report and the present one is the fact that at 
that time the numerical analysis had not been completed: in order to illustrate 
the numerical size of the two-loop contributions, Ecker had to present three 
different numerical values of the S-wave scattering lengths, depending on the 
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input used for the low energy constants (see table 1 in Ref. □). The numerical 
analysis has been recently completed and has yielded the values: 

a[; = 0.222 ± 0.005, = -0.0445 ± 0.001 . (2) 

In what follows I will describe what are the experimental ingredients and the 
theoretical tools that have lead to this result. 

2 Numerical solutions of Roy equations 

In 1971 Roy derived a set of dispersion relations for the partial wave am- 
plitudes of TTTT scattering, that fully respected crossing symmetry, and that 
involved only two subtraction constants. The crucial observation that lead to 
this result was that if one writes a fixed-i dispersion relation, crossing sym- 
metry constrains the form of the t dependent subtraction constants in such a 
way that they can be expressed in terms of the two iS-wave scattering lengths 
and other dispersive integrals. The equations have the following form: 

2 oo „oo 

where / and £ denote isospin and angular momentum, respectively and A:^(s) 
is the partial wave projection of the subtraction term, present only in S- and 
P- waves, 

H{^)=4S', + '-^^{2a",^5a^,) ^^S^, 5^ + S} - ^S^ (4) 

The kernels K^^ (s, s') are explicitly known functions. They contain a diag- 
onal, singular Cauchy kernel that generates the right hand cut in the partial 
wave amplitudes, as well as a logarithmically singular piece that accounts for 
the left hand cut. Soon after these equations became_available in the litera- 
ture, various groups started to treat them numerically a. At that time the goal 
was to derive from the available data constraints on the scattering lengths: 
all groups agreed on the result that only data sufficiently close to threshold 
could provide significant bounds on the scattering lengths. With the only 
exception of data on K^a decays, all other data were at too high an energy, 
and therefore left the S'-wave scattering lengths practically unconstrained. 

The only jhigh-statistics experiment on Kei decays provided its final re- 
sults in 1977 u, and it indeed allowed (when combined with the numerical 
studies of Roy equations) to constrain the value of the 7 = S'-wave scatter- 
ing length to within a rather narrow range: aj] = 0.26 ±0.05. Until today this 
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was the best experimental information available on this observable. Now we 
have already preliminary results 113 from a new experiment that has been oper- 
ating at Brookhavcn, analysing more than, 4 • 10^ ife4 decays (to be compared 
to the 3 • 10'* of the previous experiment □). 

The prospects of a new generation of high-statistics experiment on Kf,4, 
decays and theoretical work on the calculation of the scattering lengths made 
necessary a revival of thc_numerical treatment of Roy equations. An extensive 
work has been made in Ell, whose highlights I will now briefly summarize: 

1. Given the strong dominance of the S~ and P- waves at low energy, Roy 
equations have been solved only for these, and only on the interval 4M^ < 
s < So = (0.8 GeV)^, the lower half of their range of validity (which has 
been rigorously proven to extend up to ^si = 1.15 GeV). In that region the 
contributions generated by inelastic channels can be safely neglected. In the 
interval from sq to S2 = (2 GeV)^, the imaginary parts have been evaluated 
with the available experimental information, whereas above S2, a theoretical 
representation, based on Regge asymptotics has been used. 

2. Unitarity converts the Roy equations for the S- and P-waves into a set 
of three coupled integral equations for the corresponding phase shifts: The 
real part of the partial wave amplitudes is given by a sum of known contri- 
butions (subtraction polynomial, integrals over the region sq < s < S2 and 
driving terms) and certain integrals over their imaginary parts, extending 
from threshold to sq- Since unitarity relates the real and imaginary parts in 
a nonlinear manner, these equations are inherently nonlinear and cannot be 
solved explicitly. 

3. Several mathematical properties of such integral equations are known: In 
particular, the existence and uniqueness of the solution is guaranteed only 
if the matching point sq is taken in the region between the place where the 
P-wave phase shift goes through 90° and the energy where the / = S'-wave 
does the same. As this range is quite narrow (0.78 GeV < -y/io < 0.86 GeV), 
there is little freedom in the choice of the matching point. 

4. A second consequence of the mathematical structure of the Roy equations 
is that, for a given input and for a random choice of the two subtraction 
constants, the solution has a cusp at sq: In the vicinity of the matching point, 
the solution in general exhibits unphysical behaviour. The strength of the 
cusp is very sensitive to the value of ag. In fact, it was found that the cusp 
disappears if that value is tuned properly. Treating the imaginary parts as 
known, the requirement that the solution is free of cusps at the matching 
point determines the value of Cg as a function of Og. 

5. The input used for the imaginary parts above the matching point is subject 
to considerable uncertainties. In this framework, the values of the S- and P- 
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wave phase shifts at the matching point represent the essential parameters in 
this regard. The data on the pion form factor, obtained from the processes 
e+e~ — > TT+TT" and r — > Tr~Tr^Ur, very accurately determine the behaviour of 
the P - wave phase shift in the region of tlie p resonance, thus constraining 
the value of Sl{so) to a remarkably narrow range. The phase shifts extracted 
from the reaction ttN — > tttt A'' constrain rather strictly the difference Sq — Si, 
better than either of the two phases individually. Since the P wave is known 
very accurately from the leptonic processes mentioned above, this implies that 
5q(so) is also known rather well. The experimental information concerning 6q, 
on the other hand, is comparatively meagre. 

6. The uncertainties in the experimental input for the imaginary parts and 
those in the driving terms turn the universal curve into a band in the {uq, Oq) 
plane. Outside this "universal band" , the Roy equations do not admit phys- 
ically acceptable solutions that are consistent with what is known about the 
behaviour of the imaginary parts above the matching point. 

7. The Olsson sum rule relates the combination 2aQ — Soq of scattering lengths 
to an integral over the imaginary parts of the amplitude. Evaluating the 
integral, it was found that the sum rule is satisfied within a band that has a 
large overlap with the "universal band" mentioned above. It is by no means 
built in from the start that the two requirements can simultaneously be met - 
the fact that this is the case represents a rather thorough check of our analysis. 

8. The admissible region can be constrained further if use is made of exper- 
imental data below the matching point. At the moment there are two main 
sources of information on tttt scattering below 0.8 GeV: A few data points for 
the 1 = 2 S-wavc phase shift which will, unfortunately, not be improved 
in the foreseeable future - and a few data points on 6q — Si very close to 
threshold, from Ke4 decays. The latter provide an important constraint. 

9. The Roy equation analysis is the only method that allows one to reliably 
translate low-energy data on the scattering amplitude into values for the 
scattering lengths. As discussed above, the available data do correlate the 
value of Hq with the one of Cq. Unfortunately, however, the value of Hq is not 
strongly constrained: In agreement with earlier analyses, it was found that 
these data are consistent with any value of in the range from 0.18 to 0.3. 

10. The two subtraction constants Oq, Hq are the essential parameters at low 
energies: If these were known, the method would allow one to calculate the 
S- and P-wave phase shifts below 0.8 GeV to an amazing degree of accuracy. 
In particular one can evaluate the amplitude also at unphysical points. This 
is very important if one wants to use at best the information contained in 
the chiral representation of the tttt scattering amplitude, as we will see in the 
following section. 
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3 Matching the chiral and the dispersive representation 

The Weinberg formula for the S'-wave, 1 = scattering lengths is subject 
to substantial corrections. The physical reason for the size of these corrections 
is the strong rescattering of pions in the / = channel. The leading order 
formula fails to fully describe such effects because unitarity corrections only 
show up at next -to-leading order and beyond. On the other hand, if one uses 
the Weinberg amplitude below threshold and away from it, these unitarity 
effects will become less and less important, and the nominal parameter of the 
chiral expansion (M^ in 1 GeV^ units, a number of the order of a few percent) 
will start to dictate the actual size of the corrections. 

The fact that in the Roy equations the S'-wave scattering lengths appear 
as the subtraction constants is only a matter of choice: one could as well 
subtract at an unphysical point, without changing any of the physical results. 
A combination of the two approaches that exploits the strengths of both is 
therefore possible: the chiral expansion may provide accurate information 
on the two subraction constants, by fixing the value of the amplitude way 
below threshold, but would be less precise in the evolution of the amplitude 
up to threshold and above - there the dispersive framework based on the 
solution,of Roy equations is much better. This program has been carried out 
in Ref. Q, making full use of the information contained in the two-loop chiral 
representation of the amplitude. In the following I will briefly describe some 
of the details. 



3. 1 Chiral representation of the scattering amplitude 

The two-loop representation yields the first three terms in the low energy 
expansion of the partial waves: 

4(s)=4(s)2+4(s)4+4(s)6+0(/) . (5) 

Since inelastic reactions start showing up only at 0{p^), unitarity implies 

Im4(s) = cj{s) \t[{s)\' + 0(/) , a{s) = (1 - 4Af^/s)^ . (6) 

The condition fixes the imaginary parts of the two-loop amplitude in terms 
of the one-loop representation. At leading order, the scattering amplitude 
is linear in the Mandelstam variables, so that only the S- and P-waves are 
different from zero. Unitarity therefore implies that, up to and including 
0{p^), only these partial waves develop an imaginary part. Accordingly, the 
chiral representation of the scattering amplitude can be written as liJ 

A(s, t, u) =C(s, t, u) + 327r {i C/"(s) -I- | (s - u) U^{t) + | (s - t) U\u) 
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+ iC/2(t) + iC/2(«)-i[/2(s)} + 0(/) , (7) 
where C{s,t,u) is a crossing symmetric polynomial, 

C{S, t,u) ^ Cl + S C2 + C3 + (t — llf C4 + C5 + s (t — uf Cq . (8) 

The functions U°{s), U^{s) and U^{s) describe the "unitarity corrections" as- 
sociated with s-channel isospin / = 0, 1, 2, respectively. In view of lmtj{s)e oc 
s^, several subtractions are needed for the dispersive representation of these 
functions to converge. We subtract at s = and set 

The subtraction constants are collected in C{s,t,u). As only the S- and P- 
waves enter, we have dropped the lower index, {t'^,t^,t^} = {tQ,t\,tQ}. For 
kinematic reasons, the integrand of the P-wave differs from the one of the 
S-waves: {eo, ei, £2} = {0, ^rfi}- It is straightforward to check that the result 
of the two-loop calculation B is indeed of this structure. 

The two-loop result for the amplitude specifies how the coefhcients 
ci,..., cg depend on the quark masses cl (for explicit expressions, see ta). 
These formulae, in particular contain Weinberg's low energy theorems, which 
in this language state that the leading terms in the expansion of the first two 
coefficients are fixed by and F^: ci = -MUJ^"^ . . . , C2 = l/F^ + ... At 
first order, the constants ^1,^2,^3,^4 from £4 Ej enter, and at second order, 
the chiral representation of the scattering amplitude involves the couplings 
ri , . . . , rg from Cq Ell. It is useful to distinguish two categories of coupling 
constants: 

a. Terms that survive in the chiral limit. Four of the coupling constants that 
enter the two-loop representation of the scattering amplitude belong to 
this category: £1, £2, r^, rg. 

h. Symmetry breaking terms. The corresponding vertices are proportional 
to a power of the quark mass and involve the coupling constants i?3, ^4, 
ri, r2, r-i, r4. 

The constants of the first category show up in the momentum dependence 
of the scattering amplitude, so that these couplings may be determined phe- 
nomenologically. The symmetry breaking terms, on the other hand, specify 
the dependence of the amplitude on the quark masses. Since these cannot 
be varied experimentally, information concerning the second category of cou- 
pling constants can only be obtained from sources other than tttt scattering. 
The constants r„ from £g only generate tiny effects, so that crude theoretical 
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estimates suffice, but the couplings ^3 and £4 from £4 do play an important 
role. 

The crucial parameter that distinguishes the standard framework from 
the one proposed in Ref. t3 is £3. This coupling constant determines 
the first order correction in the Gell-Mann-Oakes-Renner-relation: = 
2Bm{l ~ \i£i + 0(C^)}, i = 2Bml{\%'K^Fl). The value of £z is not known 
accurately. Numerically, however, a significant change in the prediction for 
the scattering lengths can only arise if the crude estimate in Ref. 113, 

4 = 2.9 ±2.4, (10) 

should turn out to be entirely wrong. We do not make an attempt at reduc- 
ing the uncertainty in €3 within the standard framework, but will explicitly 
indicate the sensitivity to this coupling constant. 

Chiral symmetry implies that the coupling constant I4 also shows up in 
the expansion of the scalar radius in powers of the quark masses □: 

(0.^3;2jj{/,-i| + ?A, + 0(a} ^ (11) 

As pointed out in Ref. El, the scalar radius can be determined on the basis of 
a dispersive evaluation of the scalar form factor. The result Q 

(r^>, = 0.61 ± 0.04 fm^ , (12) 

is an update of the value given in Ref. and is consistent with earlier esti- 
mates of the low energy constant ^£4 based on the symmetry breaking seen in 
FkIF-tx or on the decay K ttIv c2l, but is considerably more accuratCpSince 
the chiral representation of the scalar form factor is known to two loops Eil, the 
dependence of the correction on the quark masses is also known. In ad- 
dition to £1, . . . ,£i, the explicit expression involves a further term, rs2, from 
Ce llZl. In the following, we use this representation to eliminate the parameter 
£4 in favour of the scalar radius. 

These numerical estimates of these low-energy constants arc important 
when analyzing the corrections to low-energy theorems. Let us look, e.g., at 
the following two relations among the coefficients ci , . . . ,04 which are dictated 
by chiral symmetry. Consider the combinations 

Ci {c2 + 4A/2(C3 - C4)} ^ C2=^ {-ci + 4M4(c3 - C4)} . (13) 

Chiral symmetry implies that, if the quark masses are turned off, both Ci 
and C2 tend to 1. The contributions from C3 and C4 ensure that the first 
order corrections only involve the symmetry breaking couplings £3 and £4. 
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Eliminating £4 in favour of the scalar radius, the low energy theorems take 
the form 

Ci = l + ^(r\ + ||+fAi + 0(a , (14) 

C2 = l + ^(r2). + ||4-i^}+C^A, + 0(a • 

At first nonleading order, Ci is fully determined by the contribution from the 
scalar radius, while C2 also contains a contribution from ^3. Inserting the 
values {r^)g ~ 0.61 fm^ and £3 = 2.9 and ignoring the two-loop corrections 
Ai, A2, we obtain Ci = 1.103, C2 = 1.117. The value of C2 differs little from 
Ci - as stated above, the estimate ( [To| ) implies that the contributions from 
£3 are very small. 

3.2 Phenomenological representation of the scattering amplitude 

If one neglects the contribution of the absorptive parts from D waves and 
higher, the mr scattering amplitude can be written as a combination of func- 
tions of a single variable that only have a right-hand cut - exactly like in the 
two-loop chiral representation. The imaginary parts that appear inside the 
dispersive integrals will be evaluated from the solutions of the Roy equations 
- since these are based on phenomenological information we will refer to this 
as the phenomenological representation. We subtract the relevant dispersion 
integrals in the same manner as for the chiral representation: 

Since all other contributions can be replaced by a polynomial, the phenomeno- 
logical amplitude takes the form 

An — 
A{s, i, u) = 16^2 + __ (2aO _ 5^2^ ^ ^ p(^^ ^^g^ 

-|-327r |iM?°(s) + |(s - u)W^{t) + |(s - t)W^{u) 

+ \W\t) + \W\u) -\w^{s)^ + 0{p^) . 

We have explicitly displayed the contributions from the subtraction constants 
flg and Uq. The term P{s,t,u) is a crossing symmetric polynomial 

P{s,t,u)=pi+p2S+p3S^ +Pi{t-u)'^ +P5S^ +pes{t-u)'^ . (17) 
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Its coefBcients can be expressed in terms of integrals over the iniaginary parts 
of the partial waves. Explicit expressions can be found in Ref. In the fol- 
lowing, the essential point is that the coefficients pi, ... , can be determined 
phenomenologically. 



3.3 Matching the two representations 

In their common domain of validity, the two representations of the scattering 
amplitude specified above agree, provided the parameters occurring therein 
are properly matched: 

47r — 
t, u) = levra^ + ^ (2a° - hal) s + P(s, t, u) + 0{p^) . (18) 

Since the main uncertainties in the coefhcients of the polynomial P{s, t, u) 
arise from their sensitivity to the scattering lengths a", Oq, the above rela- 
tions essentially determine the coefficients ci , . . . , cg in terms of these two 
observables. The same then also holds for the quantities Ci, C2 defined in 
eq. (^|). The corresponding low energy theorems for Aq and are of the 
form 

= + + Oim') , al = + Mt a, + 0{m^) , (19) 

with Co = y(12Ci — 5C2). The terms ap, a2 stand for integrals over the 
imaginary parts of the partial waves that can be worked out from the avail- 
able experimental information. Formula (^9|) clearly shows how the match- 
ing works: chiral symmetry provides accurate information on the unphysical 
quantities Cq and C2 ( p^ ) : the relation between these quantities and the scat- 
tering lengths involves the terms ao and a2, which are best calculated from 
the explicit numerical solution of the Roy equations discussed in the previous 
section. 



3.4 Results for Oq and Oq and other threshold parameters 

Inserting the one-loop prediction for Ci , C2 in the relations ( p^ ) and solving 
for Oq, Oq, we obtain the following first order results: 

a[J = 0.2195 , al = -0.0446 , 2ag - 5ag = 0.662 . (20) 

The two- loop corrections Ai and A2 involve the coupling constants 
^ij ^2, ^3, the scalar radius, as well as the terms ri, . . . , r4, rs2 from Cq. The 
size of the contributions ixppx the latter may be estimated with the resonance 
model described in Refs. d^EJ. The constants ii, (2 can then be determined 
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Roy 


CHPT 0(p«) Roy 


CHPT 0(p«) Ke4 


"o 


.220 ± .005 


.215 


.219 ± .005 


-lOag 


.444 ± .001 


.445 


.420 ± .010 




2.76 ±0.06 


2.68 


2.79 ±0.11 


-10^62 


8.03 ±0.12 


8.08 


7.56 ±0.21 




3.79 ±0.05 


3.80 


3.78 ±0.21 




5.67 ±0.13 


5.37 


5.9 ± 1.2 




1.75 ±0.03 


1.76 


2.2 ±0.4 




1.70 ±0.13 


1.72 


2.9 ± 1 



Table 1. Comparison of the numerical values of threshold parameters as obtained from 
the Roy solutions (after matching with CHPT), or CHPT. In the latter case the low-energy 
contants can be fixed either by the matching with Roy equations (third column), or by a 
phenomenological analysis of Ke4 form factors at two loops E3 (fourth column) . 



numerically with the phenomenological values of C3 and C4. The resulting 
two-loop corrections for the scattering lengths are very small. The numeri- 
cal result is sensitive to the value of the scale fi at which the renormalized 
coupling constants r^j^ip) are assumed to be saturated by the resonance con- 
tributions. In the following, we use the resonance model at the scale fj, = Mp 
and take the range 500MeV < ^ < IGeV as an estimate for the uncertainties 
to be attached to the two-loop corrections. 

A careful evaluation of the uncertainties coming from the phenomenolog- 
ical input leads to 

al = 0.220 ± 0.001 - 0.0017 + 0.027 A^2 , (21) 
al = -0.0444 ± 0.0003 - 0.0004 - 0.004 A^2 , 

with £3 = 2.9 + A4, (r^)^ = 0.61 + A,,2). Inserting the estimates ® , 

(p!^, we arrive at our final result: 

ag = 0.220 ± 0.005, -0.0444 ± 0.0010 , , . 

2ag - 5ag = 0.663 ± 0.006 , - = 0.265 ± 0.004 . ^ ' 

Having fixed the subtraction constants of the Roy equations by matching 
to the chiral representation we can now use the Roy solutions to evaluate the 
TTTT scattering amplitude at any energy. As an example we display in table ^ 
the results for other threshold parameters. For comparison we also show the 
numerical values obtained from a direct evaluation of the chiral two-loop am- 
plitude. In this case we have to specify how we fix the low-energy constants: 
since the matching procedure does provide values for the latter, this is the 
most obvious choice. On the other hand, in Ref. E3 a two-loop analysis of the 
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Ke4 form factors has lead to an alternative determination of the low-energy 
constants. In table |^ we give both sets of values for the threshold parameters. 
The comparison shows clearly that two-loop chiral perturbation theory works 
very well in describing both tttt scattering and Ke4 decays. On the other 
hand, the smallest error bars for the threshold parameters are obtained by 
combining the chiral representation and Roy equations. 



4 Conclusions 

We have reviewed the current status of the theory of tttt scattering. The 
combination of numerical solutions of Roy equations and two-loop chiral per- 
turbation theory has lead to a remarkably precise description of this reaction 
at low energy. As an example we have discussed the numerical values of the 
two S'-wave scattering lengths, which are given with an uncertainty of a few 
percent, and of other threshold parametep. with somewhat larger error bars. 

The new experiments at Brookhaven EJ and at DA<i>NE E3 on Kei decays 
will yield more precise information in the very near future. On the basis of 
the preliminary data of the E865 collaboration we have estimated in Ref . u the 
present allowed range for Og to be between 0.20 and 0.25, in perfect agreement 
with our pcediction (p^). Moreover, the pionic atom experiment under way 
at CERN t2l will allow a direct measurement of |ag — agj. These experimental 
tests will tell us whether we have reached a full, accurate understanding of 
the low-energy dynamics of pion interactions. 
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